Monitoring parton equilibration in heavy ion collisions via dilepton polarization 
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In this note we discuss how angular distribution of the dileptons produced in heavy ion collisions 
at RHIC/LHC energies can provide an information about a degree of local equilibration of the 
quark-gluon plasma produced at different invariant mass regions. 

PACS numbers: 



I. INTRODUCTION 



II. ANGULAR ANISOTROPY 
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The issue of parton equilibration in heavy ion collisions 
is an area of very active research. As it is well known, suc- 
cessful hydrodynamical description of the elliptic flows 
P~H4] implies that the beginning of (transverse) hydrody- 
namical expansion cannot start later than ~ 1/2 fm/c af- 
ter the collision moment. Perturbative mechanisms such 
as e.g. "bottom-up" equilibration discussed in [5] have 
difficulties explaining how can it happen so rapidly. On 
the other hand, applications of the AdS/CFT language 
[BHS] naturally ascribe the thermalization time to the "in- 
fall time" into an emerging black hole horizon, which is 
of the order of its position in the holographic coordinate 
~ 1/irTi ~ 0.2/m/c. 

Recent studies [El IH] have followed a set of arbitrar- 
ily chosen initial conditions through numerical solution 
of the Einstein equations. At late time a convergence 
with a hydrodynamical description is observed, as ex- 
pected. A somewhat surprising finding is that agree- 
ment with viscous hydrodynamics is reached when the 
anisotropy is still quite large. We would like therefore 
to distinguish the "hydronization" |12) time, at which 
local stress tensor T^ v agrees with hydrodynamical one 
and the "anisotropization" time, at which all distribu- 
tions become local (independent on gradients) and thus 
isotropic. (Both with a prescribed accuracy, of course.) 

This letter is not however about theory of equilibra- 
tion, but about experimental ways to monitor it in ex- 
periment. Its idea is known in general, but in this short 
note i would like to provide some numerical illustrations 
of the magnitude of the effect which can be observed in 
RHIC/LHC heavy ion experiments. 

Let us on the onset remind standard terminology to be 
used below. The sources of the dileptons are split into 
three categories: 

(i) instantaneous parton annihilation, known as the 
Drell-Yan process; 

(ii) the pre-equilibrium stage, after the nuclei pass each 
other; 

(iii) equilibrated stage, in which matter is assumed to be 
local and isotropic. 



It is well known that when spin-1/2 particles (such as 
quarks) annihilate and produce lepton pairs, the cross 
section is not isotropic but has the following form 
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(1 + cos 2 9 k ) 



(2.1) 



where the subscript correspond to a momentum k of,say, 
the positively charged lepton. This distribution is, or 
example, observed in the so called Drell-Yan pairs from 
stage (i), produced by the instantaneous annihilation of 
the quark-antiquark partons into dileptons. At high en- 
ergies the partons naturally arc collincar to the beams. 

For illustration, let us take a particularly simple one- 
parameter angular distribution 



W ~ exp[—acos P ] 



(2.2) 



with one parameter a. The subscript p reminds us that 
this angle is of the colliding partons, not final leptons. 
FigjT] shows two opposite examples of (normalized) dis- 
tributions. 

Let us now calculate the distribution of the dileptons 
corresponding to he distribution (2.2 1 



da 

dflk 



1 



-[6erf(y/a)y/na 



4er/(- v /cv)y / 7ra 
+2 v / ae _Q - erf(y/a)y/ir 
+cos 2 6 k {-Q^/ae- a + 3^erf(^) 
-2y/nerf(y/a)a)] ~ 1 + a{a)cos 2 9 k (2.3) 

The last expression is a definition of the effective parame- 
ter a(ev), which we plot in Figfll Note that large negative 
values of the a, corresponding to partons collimated near 
the beam direction and Drell-Yan process a ~ 1, as al- 
ready noticed. 

On the other hand, the second stage of the collision 
(ii) is characterized by the longitudinal pressure smaller 
than the transverse one. One may understand that be- 
cause such anisotropic parton distribution with small dif- 
ferences in longitudinal momenta is produced by a "self- 
sorting" process, in which partons with different rapidi- 
ties get spatially separated after the collision. We thus 
expect at this stage large negative a, in terms of the 
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FIG. 1: (Color online) Two examples of the angular distribu- 
tion in polar angle theta. The (black) solid curve corresponds 
to 

alpha — 10, and the (blue) dashed one to 
alpha — —10. 



parameterization we use. One then finds that the corre- 
sponding asymptotic value of the anisotropy to be 



a(a — > oo) = —1/3 



(2.4) 



III. DILEPTON INVARIANT MASS CAN BE 
USED AS A CLOCK 

Qualitatively it is clear that the most massive dilep- 
tons should come from earlier stages. The quantitative 
relation between the invariant mass and the proper time 
of the production has been studies already in the first 
paper suggesting dileptons as a diagnostic tool for QGP 
|10) . It is convenient to introduce the so called "temper- 
ature profile" function 



dN e+e - 
dM 



d x W{T{x)) = f f dTd A x5{T-T(x))W 



dTW(T)<5>{T) 



(3.1) 



where <I>(T) = J d 4 xS(T — T(x)) is the space-time volume 
in a fireball in which the temperature is between T,T + 
dT. Typically the production rate W ~ exp(—M/T) is 
exponentially decreasing at small T, while the profile, 
which can be approximately parameterized in a power 
form $(T) ~ \/T p is strongly increasing toward later 
times and smaller T. The product has a sharp maximum 
at the temperature T* = M/p 



P 



W(T)$(T)~exp[-^(T-T*) 2 ] 



(3.2) 



Its width gets small if the power p is large, which is in 
fact the case, p ~ 5 — 7 or so. The smallness of this width 
defines the accuracy of this clock. 



Finally, when equilibration is over, at stage (hi) the lo- 
cal distributions get isotropic, a — > 0, and the anisotropy 
o(0) = 0. 




FIG. 2: (Color online) Dependence of the dilepton polariza- 
tion parameter a on the parton distribution parameter a 



IV. SUMMARY AND DISCUSSION 

The Drell-Yan process (i) dominates for large mass 
dileptons above the charmonium region M > 4 GeV: 
here a is large and negative and the anisotropy parame- 
ter is a w 1. 

The preequilibrium phase (ii) mostly produce the so 
called Intermediate mass dileptons (IMD), with m^ < 
M < rnjih. here we expect large positive a and negative 
anisotropy a w —0.2. 

The later well-equilibrated stage (iii) produces smaller 
mass dileptions M < 1 GeV which tend to be unpolar- 
ized, with a — > 0, a — » 0. 

The main message of this note is directed toward ex- 
perimentalists: the dependence of the anisotropy param- 
eter on the dilepton mass a(M) turns out to be of great 
interest and thus should be measured. As we argue, it 
will change from 1 to a negative value before going to 
zero. The most negative value measured in experiment 
can be compared to our Fig(2j from which one can read 
effective angular distribution parameter a of the quark 
distribution. 

Finally, let us end with a message to theorists. Im- 
portant simplifying assumption made above is that 
anisotropy of the stress tensor can be directly translated 
into the anisotropy of parton distribution in the plasma. 



This would be true in kinetic (weak coupling) descrip- 
tion of QGP, but unfortunately strongly coupled sQGP 
is much more complex, in particular it does not allow for 
partonic quasiparticle description at all. Furthermore, 
it was emphasized [B] that in out-of-cquilibrium setting 
the average stress tensor (given by a one-point observer) 
contains information different from what is provided by 
the non-local (two-point, or spectral) observers. Basi- 
cally the same message comes out from more recent study 
[llj . It is also not obvious that spectral information in 



scalar or graviton correlators are the same as in the vector 
channel, relevant for dileptons. Thus the latter should be 
calculated. 



Acknowledgments. This note is a result of a discus- 
sion which took place at the KITP program on novel nu- 
merical methods, ADS/CFT etc in January-March 2012. 
I am grateful to KITP for the support during my stay 
there. 



[1] D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. Lett. 

86, 4783 (2001) nucl- th/0011058] . 
[2] D. Teaney, J. Lauret and E. V. Shuryak, |nucl-| 

th/0110037] 
[3] P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen 



and S. A. V oloshin, Phys. Lett. B 503, 58 (2001) hep- 

ph/0101136] . 
[4j T. Hiranoan d K. Tsuda, Phys. Rev. C 66, 054905 (2002) 

[nucl-th/0205043] . 
[5] R. Baier, A. H. Mueller, D. Scruff and D. T. Son, Phys. 

Lett. B 502, 51 (2001) |hep-ph/0009237 . 
[6] S. Lin and E. Shuryak, Phys. Rev. D 78, 125018 (2008) 

|arXiv:08 08.0910 [hep-th]]. 



[7] P. M. Chesler and L. G. Yaffe, Phys. Rev. Lett. 102, 

211601 (2009) arXiv:0812.2053 [hep-th]]. 
[8] P. M. Chesler and L. G. Yaffe, Phys. Rev. Lett. 106, 

021601 (2011) [arXiv:1011.3562| [hep-th]]. 
[9] M. P. Heller, R. A. Janik and P. Witaszczyk, 

|arXiv:1 103.3452 [hep-th]. 
[10] E. V. Shuryak. Phys. Lett. B 78. 150 (1978) [Sov. J. Nucl. 

Phys. 28, 408 (1978)] [Yad. Fiz. 28, 796 (1978)]. 

S. Caron-Huot, P. M. Chesler and D. Teaney, Phys. Rev. 



[11] 



D 84, 026012 (2011) arXiv:1102.1073 [hep-th]; 



[12] A term appeared in a discussion at KITP, to my recol- 
lection used by D. Gross first. 



